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The interaction of electrons with coherent chorus waves in the random phase approximation can be 
described as quasi-linear diffusion for waves with amplitudes below some limit. The limit is calculated 
for relativistic and non-relativistic electrons. For stronger waves, the friction force should be taken into 
account. 
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1. Introduction 

The interaction of electrons with whistler mode waves is an important process in the dynamics 
of the Earth’s outer radiation belt (RB) (7) and the RBs of other planetary magnetospheres (2, 3), 
where whistler mode waves can be found as broadband hiss (Jupiter, Saturn, Uranus, and Neptune) 
and narrowband chorus packets (Jupiter, Saturn, and Uranus) (4). The chorus packets consist of 
apparently randomly separated coherent sub-packets each with a few to few-tens of milliseconds 
duration (5). 

Two approaches are used to analyze particle energization and scattering as a consequence 
of the particle interaction with the chorus wave packets. In the first, “quasi-linear”, approach, 
particle motion is considered to be stochastic because there are multiple interactions with a large 
number of uncorrelated wave packets. It is believed that the arbitrary sub-packet sequence, medium 
inhomogeneity, and particle bouncing introduce the required randomization. The observed chorus 
wave packets are approximated by incoherent broadband noise, and the quasi-linear diffusion 
theory is applied to the global scale dynamics of the RBs (see Summers et al. (6), and references 
therein). This approach assumes small wave amplitudes, multiple interactions of the waves with 
the particles, and neglects nonlinear effects. 

In the second, “local”, approach, different nonlinear mechanisms for particle energization and 
scattering are analyzed. Particle motion in monochromatic ion cyclotron waves and whistlers (see 
Albert (7, 8) and references therein) and an inhomogeneous magnetic field is considered, and regu- 
lar as well as stochastic particle motion is found. Different regimes of particle diffusion are identi- 
fied. Faith et al. (9) showed that stochasticity of particle motion in a strong monochromatic whistler 
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wave propagating parallel to the ambient magnetic field results from particle bouncing. Roth et al. 
(10) found that bouncing along the inhomogeneous magnetic field and an oblique monochromatic 
wave are required for effective stochastic diffusion in energy and pitch-angle. Matsoukis et al. 
(77) established that an electron interacting with two oppositely directed, parallel propagating 
monochromatic whistler waves exhibits stochastic behavior and energy gain. Wykes et al. (72) 
extended the treatment to broadband whistler wave packets with continuous spectra and showed 
that the stochastic diffusion mechanism is again present. They estimated the pitch-angle diffusion 
coefficient for electrons. Khazanov et al. (13-15) considered the dynamics of non-relativistic and 
relativistic particles in a coherent whistler packet with a discrete spectrum and calculated the 
domain of particle stochastic motion as functions of particle energy and wave packet parameters, 
energy and pitch-angle diffusion coefficients, and steady-state solutions of diffusion equations. 

Cattell et al. (16) recently reported whistlers with electric field amplitudes greater than 
240 mV / m, which is larger an order of magnitude larger than any previously observed in the Earth’s 
RBs. Cattell et al. (16) suggested that the standard quasi-linear approach taken by global scale RB 
models may be inadequate to understand the dynamics associated with such large amplitude waves. 

Both conditions required for the quasi-linear description to be applicable, namely (a) small 
wave amplitudes and (b) broadband spectra with random phases, can be violated in planetary 
magnetospheres. This raises a question concerning where and when the limitations on the quasi- 
linear description of the electron interaction with whistler waves break down. 

This paper presents an analysis of the applicability of the quasi-linear approach for the inter- 
action of a strong whistler with relativistic and non-relativistic electrons based on the results 
obtained by Khazanov et al. (14). 


2. Electron-chorus interaction 


Consider a relativistic electron interacting with a circular polarized whistler wave packet 
propagating along a homogeneous ambient magnetic field. The wave packet A"’(t,z) = 
^2 k A k exp[i(fe — co k t)] will be approximated assuming that the amplitudes are constant, k = 
ko + mAk, co — cd o + tnAco, Ak/ko <-£ 1, Aco/coq 1 and the particle velocity is larger than the 
wave group velocity. This results in the space-like wave packet representation (77) 


m—N 

A"'(t, z) — A e 1( ^ exp(izm Ak) = A e ^ S ^ — — n^j , 

m=—N n 


<j) — k 0 z — cogt, 



N = 


2itko 
A k 


(1) 


Here, A w , coq, and ko are the wave packet characteristic vector potential, frequency, and wave 
vector, respectively; z is chosen along the ambient magnetic field, L is the characteristic space 
scale, and TV 1 is the number of modes in the packet. Therefore, the initial broad (N ~2> 1) wave 
packet is reduced to a monochromatic wave interacting with the particle at some points along the 
trajectory. This representation qualitatively corresponds to the observations, as can be seen from 
figure 1 of Santolik (18). 

The Hamiltonian of the problem presented in the action-angle canonical variables of the 
unperturbed (A w = 0) problem is 


H(z,p z ,0,I, t) = Ho(p z ,I) + ecy/2mco B lH 0 'A cos ^ <5 ^ — nj. 
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Hq = cjm 2 c 2 + 2 mco B I + p 
eBo 

o>b = , i js = (p — 0. 

me 


( 2 ) 


Here, tan 0 — —p x /p y and I = (p 2 + p 2 ) /2mco B . The phase 4/ is the angle between the particle 
momentum component normal to the ambient magnetic field and the wave magnetic field vector. 
The equations of motion are 


dp 

dr 

d / 
dr 

d x[r 
dr 


ecko^2mu> B lH Q *A sin i Js ^ 8 — n'j, 

n 

\]2m(ti B lH^ x A sin i/r 8 ^ — n'j. 


ec 


kop z c~ + co B mc 2 dz p z c 2 

W 0 W °’ di = ~Hq' 


with the first integral 


p z — kol = const = S. 


(3) 


(4) 


The constant S is set to zero below. The integral of motion reduces the dimension of the system 
and results after mapping into the equations of motion: 


f+i = ill 1 + Q s in Vo 


V0i+ 1 — ^ n F 1)? 


3y/2mw B 
Q = e — — AL, 


T = 


2koc 

THo,„+i 

( Pz,n+\C 2 )' 


(5) 


The phase shift F(I n+ \) is presented in Equations (3) and (4) of Khazanov et al. (14) and for 
ultra-relativistic particles reduces to F(I n+ \) = 1.5sgn/„ + i/|/„ + i|, while T = L/c. Equations (5) 
in this limit have been studied by Khazanov et al. (14), who found that these equations describe 
a strange attractor and the random walk rate can be characterized by a diffusion coefficient 


D = 


T n +i 


( 6 ) 


3/2 

with the averaging on the ensemble of particles with different initial I 0 , i>o ■ They also calculated 
the upper boundary of the particle energy as a function of wave packet parameters for stochastic 
diffusion. 

The diffusion coefficient is not sensitive to the exact expression for the phase shift. The differ- 
ence between the diffusion coefficients calculated from the general expression for the phase shift, 
F(I n+ 1 ) in Equations (5), and the simplified ultra-relativistic expression is less than 5%. 

Khazanov et al. (14) set the integral of motion (Equation (4)) equal to zero, a restriction implying 
a specific choice for the particle’s initial conditions. The choice of the integral of motion defines 
the permitted lower boundary of the momentum modulus, \p z \ , in the process interaction. 

The interaction is effective only if the resonance condition is satisfied. For a wave with a 
frequency below the cyclotron frequency, this condition holds only when the wave vector and 
particle velocity are oppositely directed. For positive k, this means negative velocities, p z < 0. 
The choice of S = 0 means that the particle can move along a straight line in the p z - I plane 
until the point (0,0) and all possible resonance velocities for the particle are available. For S > 0, 
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positive p : S are permitted, but they lie outside the resonance region. For S < 0 the region of 
negative p z S is not available to the particle. This domain is smaller than the initial p z and can 
somewhat change the diffusion time compared to the case S — 0. Therefore, both restrictions, 
the ultra-relativistic phase shift and the choice of the initial conditions leading to S — 0 do not 
affect the applicability of the results found by Khazanov et al. (14), and they should be valid for 
a relativistic particle with arbitrary initial conditions. 

To determine how well the random phase assumption of quasi-linear theory holds, the first order 
diffusion coefficient was calculated from Equations (5) and (6). We also calculated the diffusion 
coefficient in the random phase approximation from these equations. This coefficient is found by 
substituting — 1^~ from the first Equation (5) into Equation (6) and averaging over the phase 
%lr to obtain D nim i — Q 2 /2T. The diffusion coefficient found from the simulations and normalized 
by Aand is close to one for the stochasticity parameter Q in the range from 0 to 0.0025 with an 
accuracy of a few percent. Therefore, the assumption of phase randomness of the quasi-linear 
theory is valid. 

The second assumption of the quasi-linear approach is that the wave amplitude is small. This 
condition is equivalent to the assumption of particle motion along an unperturbed trajectory in the 
derivation of the quasi-linear diffusion coefficient. Corresponding to this requirement, a mapping 
and diffusion coefficient can be obtained from Equations (3), if the right-hand side in the first 
of these equations is kept constant and the random phase approximation is used. This coefficient 
can also be found from Equations (5) and (6) [Albert, private communication, 2007] by the 
perturbation method. We introduce a variable S = ( I n+ \ — l „) //„, and keep the three first terms in 
the Taylor expansion of I 3//1 in the first equation in the mapping (5) assuming that for diffusion <5 
is small. Then, this equation gives 


8 1 + t = 


3 Q 


2 /„ 


3/2 


sin i/r„. 


( 7 ) 


With the second brackets set to one the diffusion coefficient is the same as the quasi-linear diffusion 
coefficient found by Albert (79) 


Dql 


■ 


1 co B p±c B w \ 2 

2 kc H 0 B 0 ) A£|v u - v gr n| ’ 


(8) 


if the limit of relativistic particles is taken. Here, v gr is the wave group velocity. 
Substitution of S taken in the linear approximation into Equation (7) leads to 


5 1 + 


Q sin xfr,, 


61,, 


3/2 


= Q sin i/v 


3/2 

Therefore, the limit of quasi-linear approximation is defined by the condition Q/ (67„ ) <</ 1 or 


co B me 2 \ 3/2 /v ph \°-5 B w Lco 0 
co o Ho J V 2 cl Bo c 


« 1 . 


( 9 ) 


Here, v p h, B w , Bq whistler’s phase velocity, wave magnetic field intensity, and local magnetic 
field, correspondingly. This analysis can be performed for an arbitrary choice of the first integral, 
but the results remain practically the same. 

For typical parameters in the Earth’s RBs, coq/cob = 0.5, v p h = 0.25c, Hq = \ MeV, magnetic 
shell L sh = 5, and scale size L = 10 9 cm, the quasi-linear approximation is valid for B w /B q <+ 
0.005. Therefore, for whistlers comparable to these observed by Cattel et al. (76), the quasi-linear 
limit can be violated even if stochasticity develops. 
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A similar analysis for the non-relativistic particles leads to diffusion coefficient (8) for vn <-£ 
Vgrii . The restriction on the applicability of the quasi-linear approach is 


%y c b" La>Q ^ i 

wo J vx Bo c 


(10) 


For the typical Earth’s RB parameters listed above and electron energies ~15 keV, this condition 
is more severe, B"' /Bo << 1 0 4 . 

The meaning of the terms neglected in the quasi-linear limit can be clarified from the Fokker- 
Plank equation that can be written for the distribution function/fu, /, t) in the phase space region 
of particle stochastic motion. In the random phase approximation, the equation is 


df 

dt 



( 11 ) 


where u — / i/2 ; and the diffusion coefficient, calculated from Equations (5) and (6), D = Q 2 /2T, 
is independent on u. Transforming Equation ( 1 1 ) to variable /, we obtain 


a f a 

/d/\ 2 aF" 

a 

D d d A 2 


D\ — 1 — 

“I - — 

F — 

~dt ~ di 

1 

CO 

"s' 

-73 

a/ 

2 dl \dn ) 


( 12 ) 


Here,/ du = F dl and D(dl/du) 2 = /)qi. with Dq l that has been found from Equation (7) in the 
linear approximation. Therefore, the correction to the quasi-linear diffusion due to the dependence 
of the right-hand side in equation of motion (4) on the particle magnetic momentum results in the 
friction (drag) force. The role of this term depends on the characteristic scale of the problem and 
can be of the same order as the diffusion term. 


3. Conclusion 

Particle motion in the whistler wave packet given by Equation (1) can be described as diffusion. 
The diffusion coefficient can be calculated in the random phase approximation and coincides with 
the standard quasi-linear diffusion coefficient for weak waves. The limits on the wave amplitude 
that permits quasi-linear description are found. For stronger waves not only diffusion, but also 
drag force effects should be considered. This can be the case for some of the observed whistlers 
that are discussed by Cattell et al. (16). 

The results presented in this paper are based on an assumption that the observed chorus can be 
represented as a monochromatic wave and their interaction with the particles can be described as 
instantaneous kicks. The model neglects the observed frequency drift. The role of these factors 
requires additional studies. Another restriction on the applicability of these results occurs in the 
case of very strong waves when the particle dynamics is intermittent, combining chaotic, and 
regular motion (14, 15). 
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